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Abstract. 

The causal character of singularities is often studied in relation to the 
existence of naked singularities and the subsequent possible violation of the 
cosmic censorship conjecture. Generally one constructs a model in the framework 
of General Relativity described in some specific coordinates and finds an ad 
hoc procedure to analyze the character of the singularity. In this article we 
show that the causal character of the zero-areal-radius (R = 0) singularity in 
spherically symmetric models is related with some specific invariants. In this 
way, if some assumptions are satisfied, one can ascertain the causal character of 
the singularity algorithmically through the computation of these invariants and, 
therefore, independently of the coordinates used in the model. 



PACS numbers: 04.20.Gz, 04.20. Cv, 04.40.-b, 04.90. +c 
1. Introduction 

The causal character of a singularity has a well defined meaning within the theory of 
conformal boundaries [U|2}|3]. The knowledge of this causal character is fundamental 
since whenever the spacetime possesses timelike or past null singularities there are 
always null geodesies which are past incomplete. If such a singularity could develop 
from a generic gravitational collapse in the framework of General Relativity Theory, 
this would mean that the theory would lose its predictability. The question on 
whether General Relativity contains a built-in safety feature that precludes this 
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possibility was put forward by Penrose in 1969 [1] and gave rise to what is known as 
the cosmic censorship conjecture (CCC). Many counterexamples to the CCC have 
been proposed as well as many arguments in its favour (see, for example, [SJ [5] and 
references therein) so that the question of the CCC remains open. 

In this article we will deal with probably the most interesting type 
of singularities in spherically symmetric spacetimes: The zero-areal-radius 
singularities, i.e., given the areal radius R defined such that the area of a 2- 
sphere is A-kR 2 , we will be interested in singularities at R = 0. In the current 
literature the study of the causal character of the R = singularities has been 
carried out for important particular solutions. In a few simple cases the singular 
conformal boundary has been obtained by using a conformal compactification (see, 
for instance, [H1I71IE]), while in most cases there is not an analytical compactification 
and, as an alternative method, the causal character of the singularities has been 
studied through the analysis of radial null geodesies around them (see, for example, 
0, QUI dU H2])- In addition to the analysis of particular cases, this last technique 
allows some general approaches for studying zero-areal-radius singularities. In 
particular, it has led to show [13] that the causal character of the singularity 
is related to the mass function [T3]. More general studies on the formation of 
naked singularities in spherically symmetric spacetimes along these lines can be 
found in [15], [IE] and, by using ad hoc devised procedure, in [17]. Furthermore, 
in [3], by using the techniques of the qualitative behaviour of dynamic systems 
on the differential equations satisfied by the radial null geodesies we were able 
to present the most comprehensive scheme so far to try to find out their causal 
characterization taking into account, and analyzing, the possible limitations of the 
approach. However, this work was carried out in specific coordinates (the so called 
radiative coordinates) , so that its results were restricted to models described in these 
coordinates. 

Our aim in this article is to show that the causal character of the zero-areal- 
radius (R = 0) singularity in spherically symmetric models is related with some 
specific invariants. Apart from being an interesting result from a theoretical point 
of view, this coordinate independent approach means that, if some assumptions 
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are satisfied, one could find out the causal character of a model's singularity 
algorithmically through the computation of these invariants in arbitrary coordinates. 
In order to try to reach our goal we will base our approach in an analysis of the 
results in our previous article [3]. We will show that our previous results admit 
an interpretation and rewriting in terms of some invariants and we will analyze and 
explicitly state the limits for the applicability of our results coming from our specific 
approach. On the other hand, throughout the article we will use a geometrical 
approach requiring only the existence of a spacetime, but not the fulfillment of 
Einstein's equations. Thus, we just try to discover the possibilities allowed by 
this geometrical approach which includes the classical as well as the semiclassical 
framework. 

The paper has been divided as follows: In section [2] we revise well-known 
properties of spherically symmetric spacetimes and of the radial null geodesies, but 
emphasizing the corresponding degrees of differentiability for each defined object 
(what will be an important aspect for the later development of the work). In section 
[3]we analyze the relationship between general coordinate systems and the coordinate 
system used in [3]. The different cases that appear when treating the causal character 
of singularities are treated from section H] on. In particular, section 0] is devoted to 
the analysis of singularities with non-null mass function. Section [5] deals with the 
preliminaries required to the study of singularities with null mass function. Finally, 
sections El [7] and |8] analyze every null mass function subcase in detail. 

2. Basis 

Let us consider a simply connected open set V in a spherically symmetric spacetime 
and such that a part of its boundary consists of a R = interval. The metric line 
element of an oriented spherically symmetric spacetime can be (and for practical 
purposes it is usually) written in the local chart endowed with coordinates = 
x 1 , 8, (p} in the form [18] 

ds 2 = g ij dx i dx j + R 2 rffi 2 , (1) 
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where is an oriented two-dimensional Lorentzian metric (i, j—0, 1), R = R(x°, x 1 ) 
and dtt 2 = al6 2 + sin 2 6dif 2 . 

In the Lorentzian two-surface orthogonal to the 2-spheres two nonvanishing 
null vector fields may be defined such that they are linearly independent at each 
point. If some differentiability requirements are satisfied in V the integral curves 
of the two null vector fields provide us with two families (J-"i and Jfy of affinely 
parametrized null geodesies called the radial null geodesies. Take, for example, the 
geodesies belonging to the T\ family satisfying 
d 2 x a dx 13 dx 1 

~dF = ~ & d£ d£ ' 
where I is their affine parameter. The theory of Ordinary Differential Equations [19] 
together with the definition of the Christoffel symbols guarantees the existence 
and uniqueness of the affinelly parametrized geodesies provided that g a p is at least 
C 2 ~ and det(g a p) ^ 0. From now on we will guarantee the existence and uniqueness 
of affinely parametrized null geodesies by assuming that g a p is C n+1 , where n > 1— |§j, 
and that and det(g a p) ^ QJJ . 

Under these assumptions the theory of ODEs also guarantees that the solution 
of the geodesic equation will be C n+2 . Since the radial null tangent vector field to 
the T\ family I = d/di has an associated covector satisfying l[ a -p] = 0, then it can 
be written as the differential of a C n+2 function ^(x ,^ 1 ): d a u = l a . The curves 
u =constant define the trajectories of the T\ family of null geodesies in V. 

Taking into account that the scalar invariant R(x°, x l ) is a C n+l function we 
can define xi x °i xl ) = 9 a ^d a RdpR which is a scalar invariant C n function . This 
invariant is related to the C n invariant mass function (HJ EUJ EU [22] through 

m=*(l- X ). (2) 

In order to investigate the causal characterization of a R = interval we will 
consider the radial null geodesies around a point p in this interval. As we shown in 

§ As usual C l ~ means C^ -1 ' and that the (i — l)th-order derivatives are locally Lipshitz. On the 
other hand, C n+1 with n = j— means C^ +1 ^ . 

|| Note that the last requirement, the non-degeneracy of the metric, is usually assumed in the 
definition of spacetime [B] . 
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[3] this procedure requires that the radial null geodesies from at least one family, say 
J~i, reach (or leave) every point in the interval, what will be assumed in the next 
sections of this article. 

On the other hand, let us comment that provided that a R = interval 
is not reached (or left) by radial null geodesies of any family then the causal 
characterization of the interval is straightforward since this interval cannot be 
translated into a piecewise C l interval in the conformal boundary of the spacetime 
[3]. It can only be translated into a point where the boundary is not a C 1 curve 
and, thus, where there is not tangent vector properly defining its causal character. 

3. Analysis of a coordinate change 

As we mentioned in the introduction, the main goal of this article is to extend 
the coordinate dependent results presented in j3] so that they can be used 
independently of the coordinate system chosen to work with, i.e., to provide the 
causal characterization of R = O-singularities in a invariant manner. In order to do 
this, in this section we will deal with the connection between a general coordinate 
system of the type used for (CQ) and the coordinate system used in [3]. 

Lemma 3.1 Under the assumptions that the metric is C n+1 , with n > 1 — , and 
J = u,qR,i —u,iR,q^ (or, equivalently, du A dR ^ 0) there exists a coordinate 
change {x°, x 1 } — > {u, R}, where R is the areal coordinate and u is a null coordinate, 
such that the metric (T7]j can be locally written as 

ds 2 = -e^ X du 2 + 2ee 2p dudR + R 2 dQ 2 , (3) 

where e 2 = 1 and (3 = f3(u, R), x — x( u i R) are C n functions. 

This lemma is based in the fact that, since u is C n+2 and R is C n+1 , there 
is a class C n+1 map $ : (x ,^ 1 ) — > (u,R). The inverse function theorem (see, for 
example, [23]) guarantees the existence of C n+l functions f°(u, R) and / 1 (m, R) such 
that x° = f° and x 1 = f 1 provided that the jacobian determinant J is not null in 
V. Along this work we will denote the open set $(V) by U. 
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It follows that we can write the function \ as x( u i R) = xi x ° = f°{ u i R),^ 1 = 
where the chain rule guarantees that xi u iR) will & l so De & t least a C n 
function in the variables {u,R}. 

On the other hand, the condition J ^ or, equivalently, du A dR ^ implies 
that the vectors associated with these one- forms can not be parallel: l(R) ^ 03]. In 
this way, taking into account that l(R) will be C n (n > 1— ) in V we can define the 
invariant non-null constant 

£=-sign(f(i2)). (4) 

Note that it is invariant under future directed reparametrizations of I. If e — — 1 
(or +1), the expansion [6] of the null geodesies with tangent vector / is positive 
(negative, respectively) in every point of U. 

If we perform the coordinate change {x ,^ 1 } — > {u, R} then, due to the light- 
like character of the coordinate u, the metric of the spacetime ([1]) will take the 
form: 

ds 2 = -Adu 2 + IBdudR + R 2 dQ 2 , (5) 

where A = A(u,R) and B = B(u,R). The general future directed and aflinely 
parametrized null vector / and the future directed null vector k tangent to the T% 
family satisfying / • k = — 1 can be written as 

red r ( d A d\ _! 

l = ~BdR ; k= {d^ + 2BdRj C ' (6) 
where c = c(u) > depends on the affine parameter chosen for /. Clearly, 
sign(Z(i?)) = — e = — sign5, so that we can rewrite B as 

B = ee 2 ^ (7) 

with (3 = (3(u, R). If we state explicitly the relationship between the component g U R 
of the metric tensor in these new coordinates with regard to the old ones 

ee 2 ^=f\ u f, R g l3 , (8) 

where, as usual, should be understood as gij(f°(u, R), f x {u, R)), and taking into 
account that both f k (u,R) and gij{x°,x 1 ) are C n+1 functions in their respective 
variables, then the chain rule theorem implies that (3(u, R) is at least a C n function. 
% A sufficient (but not necessary) condition for l(R) ^ is \ — g a ^d a RdpR 
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On the other hand, if one evaluates the function x with this form of the metric 
and uses <^ one finds 

A = e^ X - (9) 

Thus, taking into account (J7J) and we obtain the required form for the metric ([3]) 
with the degree of differentiability stated in the lemma. QED. 

Corollary 3.1 Under the assumptions in lemma [3A\ the mass function m(u,R) = 
R(l — x(u,R))/2 is a C n function. 

Let us remark here that we do not need to find the explicit form ([3]) since we 
want to work in the original coordinates {x^}. In particular, the function f3 can be 
written with the help of ([8]), and the appropriate labeling^, as a function of x 1 }: 

J 

4. Invariants that define the causal character when lim^p m(x) ^ 

In this section we will discuss the causal character of R = in a point p = (p ,^ 1 ) 
such that lim x _ >p m(x) ^ 0, where x = (x ,^ 1 ). This condition implies that there is 
a scalar curvature singularity at p [3], so that p does not belong to the spacetime, 
but to its singular boundary. If one radial null geodesic of, say, the T\ family reaches 
p then we have the following 

Theorem 4.1 In case the spacetime metric (CP is C 2 ~ and there is a radial null 
geodesic reaching (either toward its past or its future) a R = 0- singularity at p with 
a value of its affine parameter I = l then: 

• if limj_^ m(x(l)) > 0, there is a spacelike singularity at p, 

• if liim_^ m(x(l)) < 0, there is a timelike singularity at p. 

This result can be found in our previous article [3] where it was shown in 
specific coordinates. However, taking into account that m is a scalar invariant, 
the proposition is true for other coordinate systems and it is, thus, an invariant 
result. Let us reiterate that the requirement of, at least, a C 2 ~ metric is a minimum 
+ So that "Zi'V 0. 
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assumption for the existence and uniqueness of radial null geodesies in the spacetime. 
The theorem only requires the existence of the directional limit along the radial null 
geodesic. As a corollary we have the following more applicable result: 

Corollary 4.1 In case the spacetime metric is C 2 ~ and lim x ^ p m(x) ^ 0, where p 
is such that lim x ^ p R(x) = 0, then the causal character of the R = ^-singularity at 
p is defined by the sign of the invariant mass function as the function approaches p: 

• if \im x ^ p m(x) > then there is a spacelike singularity at p, 

• if \im x ^ p m(x) < then there is a timelike singularity at p. 

5. Basis for the lim x _ 5>p m(x) = case 

This is the most involved case. It is known that if lim x _> p m(x) = the causal 
character of R = at p admits any possibility: spacelike, lightlike or timelike [3]. 
Therefore, the question is: are in this case other invariants which define the causal 
character of R = independently of the coordinate system used? In order to answer 
this, let us introduce in this section some new definitions and lemmas. For instance, 
in spherical symmetry there is an invariantly defined vector, the Kodama vector 
[24] . which is also known to possess very interesting properties (see, for instance, 
[25] and references therein): 
?_ ^9R__d_ 

where denotes the volume form associated with the two- metric gij. It satisfies 
£ 2 = ~x, so that if \ > R > 2m) the orientation of can be chosen in such 
a way that £ is a future directed timelike vector. Kodama's vector characterizes 
the spherically symmetric directions tangent to the hypersurfaces R =constant 
and provides and invariantly defined direction in which the area of the two-spheres 
remains constant [26] . For the metric (|3]) £ takes the form: 

£=e-* f 'a u , (12) 

which satisfies Z • £ < in the local chart. 
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5.1. On the extension ofm(u,R) and j3(u,R) 

We have seen that, if the assumptions in lemma 13.11 are satisfied, the functions 
m(u, R) and f3(u, R) are C n functions defined in an open set U. We will now see that 
these functions can be extended beyond the open set provided that some conditions 
are fulfilled. In particular, we are mainly interested in an extension of the functions 
around R = that, while it has not any physical meaning, will allow us to apply 
the theory of the qualitative behaviour of dynamic systems in a open set centered 
in a point of R = [3J. 

Lemma 5.1 Provided that there is a natural number N such that the limits of the 
functions m and (3 and of their ith-order derivatives, for all i < N , as every point in 
the boundary of the open set is approached, exist (and are finite) then the functions 
m and (3 admit a C N extension rh and (3. 

In order to show this it suffices to define the extended function and its derivatives 
for all points p b in the boundary of the open set U as 

QaQb _ QaQb 

d^™ {Pb) = ^d^R> m{u > R) > 

QaQb QaQb 



du a dR b ' ^ ' x^Pt, du a dR b ' 
for all integers a > 0, b > such that a + b < N. In this way the lemma can be 

considered as a simple case of Whitney's extension theorem [27] and the existence 

of C N extensions rh and /3 that coincide with m and (3 in 14 is guaranteed. We will 

denote the extended domain of definition by U. 

On the other hand, we want to work with the original coordinates {a; , a; 1 } 

and infer from here the extendibility of the functions m(u,R) and (3(u,R). We 

can do this by using the chain rule applied to the derivatives of m(x l (u, R)) and 

(3(x l (u, R)) as the boundary points are approached. For example, if one is looking 

for a C N extension, according to lemma 15.11 one needs, among others, the limit of 

m, u (u, R) in the boundary, but 

m, u (p b ) = lim m, u = lim m,» . 

Taking into account that f l , u can be written as a function of lj (= u,j) and R,j 
(J = 0, 1), the existence of the limit is guaranteed if we require that the limits of 
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li(x°,x l ), R,i (x°,x l ), and m,i (x ,^ 1 ) as every point in the boundary of the open 
set is approached exist (and are finite) and we also require that the limit of J when 
approaching the same points is not zero. 

This can be formalized and generalized for C N extensions similarly: 

Lemma 5.2 If the limits 

d a d b u , d a d b R 
hm „ s , s , , hm 



x^ Pb d(x*) a d(x l f ' x->n, d(x°) a d(x 

QaQb m QaQb £ 

lim „, „, „, , lim 



x^ Pb d(afi) a d(x 1 ) b ' x-^pt dix^dix 1 ) 1 * 
as every point in the boundary of the open set is approached exist (and are finite) 
for all integers a > 0, b > such that a + b < N and the limit of J as the same 
points are approached exists and is not zero, then the functions m(u, R) and f3(u, R) 
admit a C N extension m(u, R) and j3(u, R). 

5.2. On the differentiability of the singular conformal boundary 

In order to analyze the singular conformal boundary of a spherically symmetric 
spacetime it is possible to perform just the conformal compactification of the 
two-dimensional surface orthogonal to the 2-spheres retaining all the important 
information. This is so because, by means of a coordinate change, the induced 
Lorentzian metric or first fundamental form of the two-dimensional surface can be 
brought into a conformally flat form ds\ d = Q 2 (t, x)(—dt 2 + dx 2 ) or, equivalently, 

ds\ d = —{l 2 (u,v)dudv, (13) 

where {u,v} are lightlike coordinates (u = t + x,v = t — x). In this way, it can 
be naturally embedded in an unphysical two-dimensional Minkowskian spacetime 
(see, for instance, [2]). Despite an analytic conformal compactification can only be 
found for certain simple particular cases some information can be extracted without 
fully following the procedure. For instance, with regard to the differentiability of 
the singular boundary we have the following 

Theorem 5.1 If gij and R are C 3 functions, J / in V and the assumptions in 
lemma COI (for the existence of C 2 extensions) are satisfied and, in particular, they 
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are satisfied in a connected open interval C of the singular boundary R = in which 
m^O, then C will be C 3 in the unphysical spacetime. 

In order to show this we will study the differentiability of R = in the unphysical 
spacetime by analyzing the coordinate change which takes the metric to the form 
([TBI . To begin with, note that, under the assumptions in the theorem we can perform 
a first coordinate change from coordinates {a; , x 1 } to {u, R} such that m(u, R) and 
j3(u, R) will be C 2 functions. That this is satisfied follows from the requirement of a 
C 3 metric in the theorem and lemma [3TT1 together with corollary 13.11 Furthermore, 
requiring the fulfillment of the assumptions in lemma 15.21 implies that there will be 
C 2 extensions rh(u,R) and (3(u,R). 

We know from ([3]) that the metric of the lorentzian surface can be written in 
coordinates {u, R} as 

ds ld = R ~ (7dU > ( 14 ) 

where a = (R — 2m)du — 2eRe~ 2l3 dR. In order to rewrite this in double- null 
coordinates {u, v} we will look for an integrating factor such that 

e*cr = dv. (15) 

The integrability condition for ( !T5|) takes the form of the following first-order linear 
inhomogeneous PDE in U: 

a<S>, u +b®, R = c, (16) 

where the functions a, b and c are defined as 

a = 2Re~ 20 , 

b = e(R-2fn), 

c = AR(3, u e- 2 ^-s(l-2m, R ). 

Our assumptions imply that a and b are C 2 functions, while c is C 1 . We now 
choose the curve C (which can be described as {R = 0,u = s}, where s is 
a parameter defined in a connected open real interval) as the curve of initial 
conditions for the PDE and, in addition, we choose the initial condition on C to 
be Q(u,R = 0) = ip{u), where ip(u) is an arbitrary C 2 function. Then the tangent 
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vector to C (w = d/als) and the characteristic direction have distinct projections 
on the w,.R-plane (aw R — bw u = —2erh(u, R = 0) 7^ in C). This together with 
the fact that the PDE is nondegenerate (a 2 + b 2 7^ in an open neighborhood of 
C since fh(u, R = 0) 7^ 0) implies that the initial value problem has one and only 
one solution [28] which will be a C 2 function $(w, R). In this way (fl5|) can now be 
used to evaluate the slope of the singular boundary R(u, v) — in the unphysical 
spacetime: 

^\ R=0 = -2e^ R =Vm(u,R = 0). (17) 
du 

Since there are only C 2 functions in the right hand side of this differential equation, 
its solution v (u) will be a C 3 curve in the w,t>-plane. QED 

As a corollary, under these assumptions, C admits a tangent vector whose causal 
character determines the causal character of R — 0. In other words, if rh > the 
slope of the curve is negative so that the curve in the unphysical spacetime (therefore, 
the singularity) will be spacelike. Likewise, if rh < the singularity will be timelike. 
This just reiterates the results in theorem 14.11 which were then shown under less 
restrictive assumptions. The so far untreated case rh = in which the singularity 
can be timelike, lightlike or spacelike and in which the degree of differentiability of 
the singular boundary can be just C° under the differentiability assumptions on this 
subsection will be the subject of the rest of the article. 

6. Invariants if lim x ^ p m(x) = and lim x _> p £(m)(x) 7^ 

Theorem 6.1 In case lim a ._ 5>p m(x) = and lim x _>p£(m)(x) 7^ ; where p is a point 
in the singular boundary R = reached or left by a single null geodesic of the T\ 
family, the causal character of R = around p can be obtained through the invariants 

e, 

provided that g^ and R are C 3 ~ functions in V and that the assumptions in 
lemma I5.ifl (for the existence of C 1 extensions) are satisfied. Then the causal 
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Figure 1. Sketched characterization of the R = singularity when lrni x ^ p m(x) = 
and lim^-jp £(m) (x) = X\ ^ 0, where the R = 0-singularity at p is chosen 
to be reached by the radial null geodesic u = 0. As explained in the text, 
A = (1 - 2jg) 2 - Well, T = e(l — 2J°) and 8A = -teL\. In this case we 
have mJ^ = o 7^ for m / 0. In this way, if the requirements in theorem 15.11 arc 
satisfied the singular boundary should be C 3 for u^O. (In these sketches we just 
draw straight lines for u =/= instead of curves). However, in u = 0, where m = 
the singular boundary can be just C . Consider, for example, the first sketch in 
the e = — 1 column where a R = 0- spacelike singularity for u < must be abruptly 
followed by a lightlike singularity for u = 0. 



characterization around p is inferred from e, 

and Ji 



X\ = limX 1 

x— >p 



°-\imJ° 



according to figure U\ 



Note that the invariant J° is independent of the parametrization of I. 

In [3] (sec. 5) the causal characterization in this case was obtained in radiative 
coordinates provided that m(u, R) and (3(u, R) were C 2 ~ functions. That this 
requirement is satisfied follows from the requirement of a C 3 ~ metric in the 
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theorem and lemma 13.11 together with corollary 13.11 Another requirement was 
that there should be C 1 extensions m(u,R) and (3(u,R), hence our requirement 
that the assumptions on lemma 15.21 for the existence of C 1 extensions should be 
satisfied. On the other hand, it was shown that, in radiative coordinates, the causal 
characterization depends on_tl 

A = (1 - 2m, R ( U = 0,R= 0)) 2 - lQerh ru (u — 0,R— o) e - 2/5(u=0 - R=0) , 
5 A = -4em, u ( U = 0,R = o) e - 2 ^=°> B =°), 
T = e(l — 2rh, R (u — 0,R— 0)), 

where we have chosen to analyze the R = 0-singularity around u = which describes 
the null geodesic reaching or leaving p (what is always allowed through a redefinition 
of u = u zd+contant). 

— * 

On the other hand, using the expression ([6]) for I and the Kodama vector 
(jl~2l it is easy to verify that, in radiative coordinates, the invariant X 1 can be 
written as e _2/3 m m while the invariant J Q is simply m, R . Therefore, taking the 
limit {u — y 0,R — y 0), we can rewrite the quantities A, 5 A and T in an explicit 
invariant form as 

A = (1 - 2J b ) 2 - IQell 
6 A = -AeXl, 
T — e(l — 2J b °) 

Finally, in j3] it was shown that the causal character can be read from a table 
(figure [[]) depending on these quantities. In this way, one first computes m, £ and 
/ in arbitrary coordinates in order to get from them the invariants e, X\ and J b . 
Then, if the assumptions in the theorem are satisfied, these invariants provide us 
with the causal characterization of the singularity for this case. 

f Note that in [3] we chose (3(u, 0) = for easiness through the use of a coordinate change 
u 1 = u'(u). Wc do not perform this coordinate change here, what causes "exp(— 2/3(u = 0, R = 0))" 
factors in this article when compared with the corresponding expressions in [3]. 
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7. Invariants if \im x ^ p m(x) = and lim x _>p £(ra) (x) = 

— * — # — * 

Theorem 7.1 In case \im x _> p m(x) = 0, linx r _ >p £(m) = and \\m x ^ p l{m) / 1{R) ^ 
1/2, where p is a point in the singular boundary R = reached or left by a single null 
geodesic of the T\ family, the causal character of R = is given by the invariants e, 

n times 



T = CM =Cd a (- ■ ■ {Cd w (m))-) and J° 



l(m) 



l(R) 

provided that g^ and R are C n+1 functions in V (n > 2) and that the assumptions 
in lemma \5T^ (for the existence of C n extensions) are satisfied. Then e, 

1£ = limX" and j£ = lim J°, 

where ^ Q whereas X£ = for < k < n (i.e., we assume that there is a finite 
n such that it is the lowest value satisfying I£ ^ ), provide us with the causal 
characterization of the singular boundary around p according to figure 

In [3] (sec. 6) the causal characterization in this case was obtained in radiative 
coordinates provided that m(u, R) and /3(u, R) were C n (n > 2) functions admitting 
C n extensions. That this is satisfied follows, as in the previous theorem 16.11 first, 
from the requirement of a C n+l metric in this theorem and lemma [37T1 together with 
corollary 13.11 and, second, from the requirement that the assumptions in lemma 15.21 
for the existence of C n extensions fh(u,R) and j3(u,R) are satisfied. On the other 
hand, in [3] it was shown that, in radiative coordinates, the causal characterization 
in this case depends orj§ 

(O n fn \ 
e -g-^( u = 0,R=0)\ and (18) 

dm . . , 

_(« = 0,R = 0), (19) 

where d k rh/du k (u = 0, R = 0) = for < k < n (n > 2). But again 
these quantities are invariant. The second quantity ( |T9l) has already been shown to 
X Note that the definition of the case implies I® = T\ = 0. 

§ Here it has been taken into account that, when compared with the results in [3], the calculation 
using (3(u, R = 0) ^ provides just an extra factor exp(— 2/3(0, 0)) > to A ra which does not affect 

EH). 
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Figure 2. Characterization of the singularity when Ymv x ^ p m(x) = 0, 
rim x _^p£(m) = and lim x ^ p l(m)/l(R) ^ 1/2. 
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correspond with J®. With regard to the first one, let us consider the simplest case 
Xl ^ 0. In radiative coordinates the invariant X 2 is 

Z a = e-^[(e-^) )U m >tt +e-^m >uu ]. 
Since X\ = e- 2 ^°' 0) rh m (0, 0) = 0, in the limit as (u -> 0, R -> 0) we will have 

X 6 2 = e- 4 ^°)m, un (0,0). 
Likewise, in case 2^ = for < k < n and Xg ^ 0, we would have for n > 2 

-7-n _ -2n/3(0,0) d n m m 

Therefore we can write the two quantities ffl8|) and ffl9|) in an explicit invariant form 

sign(A n ) = sign(eX") and jg. 

Since these quantities define the causal character of R = in this case according to 
figure [2] [3] , then now the causal character has been determined invariantly. 

Theorem 7.2 In case Hm x -+ P m(x) = 0, lim x ^ p ^(m) = and lim x ^ p l(m)/l(R) = 
1/2, where p is a point in the singular boundary R = reached or left by a single null 
geodesic of the T\ family, the causal character of R = is given by the invariants e, 

V = f(m) and J k = £* [ 

provided that gij and R are C n+1 functions and that the assumptions in lemma \5. 6 A 
(for the existence of C n extensions) are satisfied. Then the causal characterization 
can be found according to figures^ and\Q through the computation of 

XI ee limr and J b k = lim J k , 

x—>p x—tp 

where X\ ^ whereas X 3 h = for < j < i, J k ^ whereas = for 1 < j < k 
and n =Max{i, k + 1} is supposed to exist ant to be finite. 

In [3] (sec. 6) the causal characterization in this case was obtained in radiative 
coordinates provided that m(u, R) and /3(u, R) were C n (n > 2, n =Max{z, k + 1}, 
where i, k are defined in fj20l) and (I2"T1) . respectively) functions admitting C" 1 
extensions. That this is satisfied follows, first, from the requirement of a C n+1 
metric in this theorem and lemma 13.11 together with corollary 13.11 and, second, 
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Figure 3. Characterization of the singularity when lim^-^p m(x) = 0, 
lim^-^p ^(m) = and lim x ^ p l(m)/l(R) = 1/2 and i is odd. In the extra- conditions 
we use p = (i — l)/2 and A from (|22|) . 
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Extra-conditions 



£=-1 



£= + 1 



k>p 



/V\\VAAAVVWA<A 

R=0 



R=0 

<$£n, 
w\ V. V\ 



k>p 



k < p 




R=0 




k < p 



R=0 



oi. 



k < p 




R=0 




k < p 



R=0 



Figure 4. Characterization of the singularity when \im x ^ p m(x) = 0, 
lim^^p^m) = and \im x ^. p l(m)/l(R) = 1/2 and i is even. In the extra- 
conditions we use p = i/2. 
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from the requirement that the assumptions in lemma 15.21 for the existence of 

C n extensions fh(u, R) and /3(u, R) are satisfied. According to [3] the causal 

characterization in this case and for radiative coordinates depends on sign(aj), 

sign(6 fc ) and A = b\ + 2(i + l)aj, wherq|| 

a . = _ £e -2^(o,o) 2 t + 1 &fh = = ^ 

i\ du l( J ( J 

b k = - ee -W,o>^?^ (u = 0, R = 0), (21) 

with d j m/du j (u = 0, R = 0) = for < j < i and = 0, R = 0) = for 

1 < j < k. But again these quantities are invariant. The first quantity sign(a,) has 
already been treated in the previous theorem where we showed that it corresponds 
to sign(— ell). With regard to the second one sign(6fc), let us consider the simplest 
case Jl ^ 0. In radiative coordinates the invariant J x can be written as 

J 1 = e~ W va, Ru 
which in the limit as (u — > 0, R — > 0) provide us with 

^ = 6-^^(0,0), 
which is clearly related to b\ fl2TT) . If = we should consider 

J 2 = e- 2p [(e- w ), u m, Ru +e- 2p m }Ruu ], 
which in the limit as (u — > 0, R — > 0) provide us with 

J 6 2 = e-* o 'm lta (0,0). 
Likewise, for k > 2 we will have 

Jt = e-^^M 



to be compared with (I2T1) . In this way we can write the two quantities ( 1201) and ( IzTj ) 
and, thus, their signs in an explicit invariant form 

sign(aj) = sign(-eX^) and sign(6 fe ) = sign(-ej b fc ). 

|| Note again that, unlike in [3], factors exp(— 2^(0, 0)) appear since in this work we do not choose 
/3(u, 0) = 0. It has to be taken into account that if one now wishes to obtain at and bk as in [3] 
there is a corresponding slight modification in the change of variables used there for this nilpotent 
case: x =exp(-2/?(0, 0))u/2. 
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Finally, taking into account its definition, A can also be written in an explicit 
invariant form 



Since these quantities define the causal character of R = in this case according to 
figures [3] and HJ then now the causal character has been determined invariantly for 
this case. 

8. Invariants if lim^p m(x) = Vp G C 

Theorem 8.1 If lim x _+ p m(x) = Vp in a connected open interval C of R = then 
the causal character of p is determined by the invariant 



exists (and is finite) for all p G C. Then the causal characterization around p is 
inferred from 

• If J® < 1/2 then R = is timelike in this interval. (This case includes both 
the possibility of a singular or a regular interval C). 

• If Jb > 1/2 then R = is spacelike in this interval. 

• If Jib = 1/2 then R = is lightlike in this interval. 

In [3] (sec. 7) the causal characterization in this case was obtained in radiative 
coordinates provided that m(u,R) and (3(u,R) were C 2 ~ functions admitting C 1 
extensions fh(u, R) and j3(u, R). That these requirements are satisfied thanks to the 
assumptions in the theorem has already been shown for theorem 16 . 1 1 On the other 
hand, in [3] it was shown that, in radiative coordinates, the causal characterization 
depends on m,# that we have shown can be written in an explicit invariant form as 
J° . In this way, the relationship between its value and the causal characterization 
follows directly from the results in [3] (sec. 7). 




(22) 




admit C l extensions according to lemma 15.11 and that the limit 



J b ° = lim J 1 



o 
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In this article we have shown that, provided some assumptions are satisfied, the 
causal character of the R = O-singularities in spherically symmetric spacetimes 
depends on some specific invariants. This allows to deduce the causal character of 
the singularity algorithmically. Basically, one starts with the knowledge of the areal 
radius R, the mass function m, a tangent vector field to the radial null geodesies I 
and the Kodama vector field £. From here one should compute the invariants 



where i,k — 0,1,2, ... and the exact last value to be computed is determined by the 
values of the lowest order invariants as the singularity is approached {X-l and J b , 
respectively) in the following manner: 

• If 2^ = fhb 7^ this value suffices to characterize the singularity according to 
theorem 14.11 provided the spacetime is at least C 2 ~ (i.e., if the existence and 
uniqueness of radial null geodesies is guaranteed). 

• If 2^ = nib = in an isolated point p in R = then different cases appear: 

— If X\ ^ then we will also need J b . The causal characterization around the 
singular point can be inferred from figure (TJ if the assumptions in theorem 
16.11 are satisfied. 

— If X\ = and J b ^ 1/2 then one needs to compute X\, where we demand 
the existence of a finite % > 2 such that it is the lowest value satisfying 
X l b 7^ 0. The causal characterization around the singular point can then be 
inferred from figure [2} if the assumptions in theorem 17.11 are satisfied. 

— If X\ = and J b = 1/2 then one needs to compute X\, where we 
demand the existence of a finite i > 2 such that it is the lowest value 
satisfying X l b ^ 0, and to compute J h k , where we demand the existence of 
a finite k such that it is the lowest value satisfying J h ^ 0. The causal 
characterization around the singular point can be inferred from figures [3] 
and HI if the assumptions in theorem 17.21 are satisfied. 




and 
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• If X° = mft = in an open interval of R = then the causal characterization of 
the singular interval can be deduced from sign(j7ft° — 1/2) according to theorem 
18.1^ if the assumptions in the theorem are satisfied. 

Note that for every case some assumptions must be satisfied. These assumptions 
come mainly from the fact that our results are based on the use, in j3], of the 
qualitative theory of dynamic systems to the differential equations satisfied by the 
radial null geodesies. The application of the appropriate theorems to the analysis of 
these differential equations requires some degree of differentiability for the functions 
m and (3. (More details on this issue can be found in [3]). Likewise, the reader can 
consult [3] for some applications of this technique to the study of the generation of 
naked singularities or black hole evaporation. 

Acknowledgements 

We would like to thank J.M.M. Senovilla and Conan Wu for helpful discussions. We 
would also like to acknowledge the Generalitat de Catalunya (grant 2009SGR-00417) 
for financial support. 

References 

[1] Penrose R 1963 Phys. Rev. Lett. 10 66 

[2] Garcia-Parrado A and Senovilla J M M 2005 Class. Quantum Grav. 22 Rl 
[3] Fayos F and Torres R 2011 Class. Quantum Grav. 28 215023 
[4] Penrose R 1969 Riv. Nuovo Cimento 1 252 
[5] Penrose R 1999 J. Astrophys. Astr. 20 233 

[6] Hawking S W and Ellis G F R 1973 The large scale structure of space-time. Cambridge: 
Cambridge University Press 

[7] Volovich I V, Zagrebnov V A and Frolov V P 1976 Teoret. Mat. Fiz. 29 191 

[8] Hiscock W A, Williams L G and Eardley D M 1982 Phys. Rev. D 26 751 

[9] Kuroda Y 1984 Prog. Theor. Phys. 72 63 
[10] Eardley D M and Smarr L 1979 Phys. Rev. D 19 2239 
[11] Christodoulou D 1984 Commun. Math. Phys. 93 171 
[12] Ori A and Piran T 1990 Phys. Rev. D 42 1068 
[13] Hayward S A 1996 Phys. Rev. D 53 1938 
[14] Misner C W and Sharp D H 1964 Phys. Rev. B 136 571 



On the invariant causal characterization 



21 



[15] Lake K 1992 Phys. Rev. Lett. 68 3129 

[16] Singh T P 1999 Class. Quantum Grav. 16 3307 

[17] Giambo R, Giannoni F, Magli G and Piccionc P 2003 Class. Quantum Grav. 20 L75 

[18] Plcbariski J and Krasiriski A 2006 An introduction to General Relativity and Cosmology 

Cambridge: Cambridge University Press 
[19] Arnold V I 1992 Ordinary Differential Equations Berlin Heidelberg: Springer- Verlag 
[20] Hernandez W C and Misner C W 1966 Astrophys. J. 143 452 
[21] Cahill M E and McVittie G C 1970 J. Math. Phys. 11 1382 
[22] Zannias T 1990 Phys. Rev. D41 3252 

[23] Boothby W M 1986 An introduction to differentiable manifolds and Riemannian geometry 

Florida: Academic Press Inc. 
[24] Kodama H 1980 Prog. Theor. Phys. 63 1217 
[25] Abrcu G and Visser M 2010 Phys. Rev. D 82 044027 
[26] Bcngtsson I and Scnovilla JMM 2011 Phys. Rev. D 83 044012 
[27] Whitney H 1934 Trans. Amer. Math. Soc. 36 63 

[28] Courant R and Hilbcrt D 1989 Methods of Mathematical Physics New York: John Wiley & 
Sons 



